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IDENTIFICATION  OF  SMALL  INHOMOGENEITIES  OF  EXTREME 
CONDUCTIVITY  BY  BOUNDARY  MEASUREMENTS: 

A  CONTINUOUS  DEPENDENCE  RESULT 


AVNER  FRIEDMANf  and  MICHAEL  VOGELIUSt 


Abstract.  We  consider  an  oloctrost at ic  problem  for  a  conductor  consisting  of  finitely  many  small 
inliomogeneities  of  extreme  conductivity,  embedded  in  a  spatially  varying  reference  medium.  Firstly  we 
establish  an  asymptotic  formula  for  the  voltage  potential  in  terms  of  the  reference  voltage  potential,  the 
location  of  the  inliomogeneities  and  their  geometry  .  Secondly  we  use  this  representation  formula  to  prove 
a  Lipschitz  continuous  dependence  estimate  for  the  corresponding  inverse  problem.  This  estimate  bounds 
t  he  difference  in  the  location  ami  in  the  relative  size  of  two  sets  of  inhomegeneities  by  the  difference  in  the 
boundary  voltage  potentials  corresponding  to  a  fixed  current  distribution. 


Ijl  Introduction  and  statement  of  the  main  result.  The  determination  of  con¬ 
ductivity  profiles  from  knowledge  of  boundary  currents  and  voltages  has  recently  received 
;i  lot.  of  attention  in  the  literature.  In  the  biomedical  community  a  common  name  appears 
to  have  emerged  for  such  work:  electrical  impedance  imaging  [6].  It  is  usually  assumed 
that,  the  (direct  current)  voltage  potential  u  satisfies  the  differential  equation 

V  •  (q(x)Vu)  =  0  in  R,  f7cR",n>2 

where  ~)(x)  is  the  positive,  real  valued  conductivity,  to  be  determined.  The  extra  informa¬ 
tion,  based  upon  which  it  is  sought  to  determine  7(x),  consists  of  knowledge  of  currents 
On 

■)—  and  the  corresponding  voltage  potentials  u  at  the  boundarv,  <9fL 
uv 

Let  A-,  denote  the  linear  operator  from  Hxl2(0Sl)  into  H~l/2( dfl)  which  takes  Dirichlet- 
t,o  Neumann-data: 


du 

A -.(</>)  —  7  —  with  V-(t Vu)  —  0  in  Q,  u  —  $  on  OQ  . 
au 

Complete  knowledge  of  \y  is  known  to  determine  the  function  7  uniquely  under  quite 
general  assumptions:  This  was  verified  for  analytic  and  piecewise  analytic  7  in  [11]  and 
[12],  for  7  and  dimension  n  >  3  in  [IS],  and  under  the  assumption  that  7  be  sufficiently 
close  to  a  constant,  it  was  verified  for  C°°  7  and  dimension  n  =  2  in  [17]  For  a  layered 
conductoi ,  i.e.,  7  =  7(.r  1 ),  two  sets  of  Dirichlet-  and  Neumann-data  suffice  to  determine 
even  a  bounded  measurable  7  ( of.  [13]). 
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Because  of  inevitable  error  in  measurements,  it  is  at  least  as  important  to  study  the 
continuity  of  the  dependence  of  7  on  A7  as  it  is  to  verify  uniqueness,  Alessandrini  [2] 
has  recently  examined  this  question  for  11  >  3  and  shown  that  the  mapping  “»  7  has 
a  modulus  of  continuity  of  logarithmic  type,  provided  7  is  a  priori  known  to  belong  to 
a  bounded  set  in  some  Sobolev  space;  his  proof  exploits  the  continuity  of  the  mappings 
A-,  — ♦  7|,-,st  and  A7  — >  d')/Or |;«j  established  in  [19],  The  assertion  of  Alessandrini  is,  more 
specifically,  that  if  A7l  and  A72  deviate  by  6  (in  the  operator  norm  on  D( H1/2,  H~ 1  )) 

then  7!  and  72  deviate  at  most  by  C  (log  7)”*  (in  for  some  0  <  a  <  1.  While  such 

o 

a  result,  is  theoretically  very  interesting,  it  is  at  the  same  time  somewhat  disappointing, 
since  it  predicts  quite  a  weak  form  of  continuous  dependence.  It  does  not  explain  the 
apparent  practical  success  of  various  numerical  algorithms  to  recover  7  from  only  partial 
knowledge  of  A7  (cf.  [4],  [14],  [20]  and  [21]).  To  bridge  the  gap  it  seems  relevant  to  analyze 
the  continuous  dependence  for  interesting  classes  of  conductivities,  where  the  functional 
dependence  on  r  is  further  restricted;  one  example  of  such  analysis  is  found  in  [5],  [9],  It  is 
of  practical  importance  to  seek  continuous  dependence  estimates  in  terms  of  only  finitely 
many  sets  of  Dirichlct-  and  Neumann-  data. 

In  this  paper  we  consider  conductivities  that  correspond  to  a  finite  number  of  small  in- 
homogeneities,  with  extreme  conductivity,  imbeded  in  an  n— dimensional  reference  medium. 
The  reference  conductivity  ~f(x)  satisfies 

(1.1)  0  <  c0  <  7(z)  <  C0  <  00. 

We  assume  that  each  inhomogeneity  has  the  form  z jt  +  tp^B  where  B  is  some  bounded 
domain  in  R’1  with 

1 1.2)  0  G  B  and  OB  of  type  C2+fi  ,  for  some  0  <  0  <  1 

The  points  belong  to  and  satisfy: 

i  \sk  ~  Zj  I  >  d0  >  0  ,  Vj  ^  k  ,  and 

dist .(zk,dQ)  >  d0  >  0  ,  Vfr. 

The  parameter  f  determines  the  common  length  scale  of  the  inhomogeneities,  and  the  pk  . 
(1.4)  d0  <  Pk  <  D0  , 

determine  their  relative  size.  We  always  assume  that  e  is  small  enough  that  the  sets 

cl  0 

: k  1  2 fpkB  are  disjoint  and  that  their  distance  to  R"\fi  is  larger  than  —  . 


=  [J(ca-  +  epk  B) 
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denote  tin.*  total  collection  of  inhomogeneities.  If  they  all  have  infinite  conductivity,  then 
the  voltage  potential  v(  =  u’f’  ,  given  the  boundary  current  t/\  is  the  solution  to 


nun 
ti€  // 1  (S2 ),  Vh 


in  <  —  [  1  \  Vu  |2  dx  —  f  il'vd. s  » 

r«=o  in  u,  I  2  y  J 

\  a  Mi 


If  all  the  inhomogeneities  have  conductivity  0,  then  the  voltage  potential  u(  =  u“  solves 


min 

U£imn\ 


Jt  / 


|  Vi/ 1 2  e/x  —  /  i/u«/s 


Wo  assume  that 


12  is  bounded  with  <9f2  6  C2+/*  . 


R  KM  A  UK  1.1.  For  problems  (1.5),  (1.0)  to  have  a  solution  it  is  necessary  and  sufficient 
that  //■  G  H~'/2{dn)  with 

(1.8)  fipds  =  0. 


We  shall  make  the  solution  u  =  of  (1.5),  or  (1.6),  unique  by  requiring  that 


J  U'ds  =  0  . 


The  solutions  and  may  be  obtained  as  limits  from  an  electrostatic  problem  involving 
finite  and  nonzero  conductivity:  Let  7Q  -  7ln\w,  +  >  0,  and  denote  by  the 

solution  of 


mm 
«£//’(  H) 


<  \  dx  -  j  4'uds  > 

,  u  dn  ) 


t* *>Y 

^SPECTED 


Then 


(1.11) 


h®  =  lim  and  u°  =  lim  u "  , 

‘  0  —  00  Or— 0 


the  first  limit  is  relative  to  H  ’(!)),  the  second  relative  to  /f1(Q\u.’f).  _ 

In  most  of  this  paper  we  shall  focus  our  attention  on  the  problem  with  inhomogeneities  - . — 

of  infinite  conductivity;  the  case  of  zero  conductivity  is  very  similar  and  is  treated  briefly  '  - 

in  the  final  section.  --•  —  ** ■/  Codes^ 

;  |Avall  and/or 

o  Dl3t  !  Special 


\Ye  shall  henceforth  write'  u,  instead  of  ii?°.  Our  main  tool  of  analysis  is  to  express  a, 
in  terms  of  V ,  tin'  solution  of 


(1.12) 


r 


min 

ii 


Vk 


r  is  the  voltage  potential  corresponding  to  the  reference  medium  alone,  and  it  is  normalized 
I  >y 


(1.13) 


=  0 


\Ye  assume  that 

(1.14)  76  C2+^(fi), 

and  that  tin-  following  non- degenr  racy  condition  holds: 

(1.1."))  VU(x)^  0  V.r  6  Q  • 

In  practice  this  means  that  we  cjui  choose  any  y-harmonic  function  U.  with  VU(x)  ^  0  in 
il,  and  then  apply  the  houndary  current  V’  —  ~)0U /Ou  on  dQ. 

Consider  two  arbitrary  collections  of  inhomogeneities 

l<  A" 

-  =  U'-i  +  epk  D)  and  lo\  =  (J(4  +  «/>*#) 

k=\  *=i 

hotli  satisfying  (1.3),  (1.4)  and  denote  by  u(  and  u'  the  corresponding  voltage  potentials 
(with  fixed  boundary  current,  (/•)•  Our  main  result  is  the  following  continuous  dependence 
t  heorem: 

Til  non  CM  1.1.  Let  r  be  a  given  nonempty  open  subset  of  DQ.  There  exist  constants 
0  i\,  and  C  and  a  function  i /(e)  .  bin  t/(e)  =  0,  such  that  if  e  <  e0  and 

f  —0 

<  ■"  "  ||(/,  -  <  <\>  then 

(i)  K  —  K',  and.  after  appropriate  reordering , 

(H)  I  - A-  —  -vl  -T  \i>k~n'k\ 

<  C  f  "||a,  -  ||  /,  -  ( |)  +  i/(f)  (1  <  k  <  h  )  . 
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The  constants  c (,,<S„  and  C'  and  the  function  1/  depend  on  do,  Du  T.  ft,  7  and  ij>  hut  are 
otherwise  independent  of  the  two  sets  of  inhomogeneitics . 


The  factor  f~n  in  front  of  ||«t  —  t/'  ||  (f)  is  best  possible;  it  follows  immediately  from 

our  analysis  (cf.  Lemma  3.3)  that  even  if  |ct  —  z[.\  and  | pk  —  p'k\  is  of  order  1  then  the 
discrepancy  in  the  boundary  data  is  of  order  en.  For  fixed  (and  small)  e  our  theorem 
basically  shows  that  the  locations  of  the  inhomogeneitics,  r*,  and  their  relative  sizes,  pr, 
depend  Lipschitz-continuouslv  on  the  rescaled  boundary  deviation  e~n  ||uf  —  t^||i»(rp 

In  the  formulation  of  the  theorem  we  used  the  L°°  norm  of  {u(  —  i/')|r  ;  this  is  not 
essential,  in  fact  the  analysis  in  Section  4  shows  that  other  norms  can  be  used,  such  as  the 
L]  norm. 

A  brief  outline  of  this  paper  is  as  follows.  Theorem  1.1  is  proved  in  section  4.  The 
proof  is  based  on  an  asymptotic  representation  formula  for  u( ,  as  derived  in  section  3.  To 
establish  this  representation  formula  we  require  some  energy  estimates  of  U  —  ?/,;  these 
estimates  are  found  in  section  2. 

R  KM  ARK  1.2.  Recently  Friedman  [S]  has  shown  that  the  presence  of  an  inhomogeneity 
or  a  collection  of  inhomogeneities  can  be  detected  by  boundary  measurement  of  the  voltage 
potential  (on  T)  corresponding  to  a  single  current  distribution.  For  the  case  of  small 
inhomogeneitics  of  extreme  conductivity  Theorem  1.1  goes  much  further.  It  shows  that 
the  exact  locations  and  the  relative  sizes  arc  determined,  in  a  continuous  fashion,  by  the 
same  measurement. 

R  IlMAHK  1.3.  Theorem  1.1  can  be  extended  to  the  case  where  each  inhomogeneitv  is 
of  the  form  Zk  -fi  epkBk  and  D k  are  different,  fixed  domains  with  smooth  boundaries. 

it'2  Energy  estimates.  In  this  section  we  estimate  the  difference  between  ?/,  and  the 
reference  potential  U.  The  first,  result  concerns  the  i/1(fI)-norm  of  V  —  v(. 

Ll  MMA  2.1.  There  exists  a  constant  C'  such  that 


i 


3 


I ii  w-u<)  12  +  i  u-uf 
h 

Proof.  Since 

I  |  r  |J  dx  <  C  {  I  | 
u  u 

and  f  (l  —  ii,  )  d.s  =  0  .  it  suffices  to  prove  t 
au 


I2)  dx  <  C 


1- 


ds  j-) 


an 


*  .1 


1 


(2.1) 


I  |  V(V  —  ut)  |2  dx<Cen\ml-l/a(on). 
n 

V,  =  :  Vi>  =  0  ill  U>„  f  vds  =  0}  . 


From  the  definitions  of  u(,U  we  get 


(2.2)  J  7Vu,  •  Vui  dx  —  J  V'U'  ds  =  j  VC  •  Vie  dx  Vie  €  Vf  , 

u  an  n 

wliicli  implies  1  liat  v,  is  the  projection  of  U  onto  the  space  Therefore 

(2.3)  I  7  |  V(C  -  «<)  |2  dx  =  min  j  7  |  V(C  -  v)  |2  dx  . 

tt  '  n 

Since  7  is  hounded  from  above  and  away  from  0,  it  suffices  to  show  that  there  exists  a 
e,  £  l  such  t  hat. 


I  I  V(C  -r<)  I2  <r/.r  <  C fn||F 


//->/J(an)  • 


Consider  first,  the  case  of  one  inhomogeneity  (with  r  =  0  and  p  —  1).  Let 


7  =  - - -  [  Udx  . 

(2e )"  |  D  |  J 


A  rescaling  of  the  Poincare  inequality  (on  2B)  gives  that 


V  -  c(  |2  dx  <C  r  /  |  VC  |  dx 


Define 


ir,{x)  —  o(-)  c,  -  ( l-o(-))l  , 
(  f 


.■here  o  is  a  C 1  cutoff  function: 


1  fi  >r  ij  €  D 

0  for  1/  E  R"\2  D  . 


Clearlv 


r,  ii 


I  M  I 


I r  If,  (Is 


is  ;iu  element  ol  From  (2.C>),  (2.7)  it  follows  1 1 1 : 1 1 


J  I  V(t-  -r,)|2  il.r  =  J  I  V(V  -’(•,)  I*  d.r 
12  $2 

<  J  I  V£’  I2  ,l.r  +  J  I  V(U  -  tr, )  I2  r/.r 


2<  ii\<  n 


and  thus,  using  (2.5),  (2.G).  we  got 


f  I  V(r-r,)|2  d.r  <C  I  I  VL  |2  fl.r 


Bv  ellipt  ic  estimates 


j  |  VA  |2  d.r  <|  Z?  |  2V'||VA||,~(,f/<)  <  C*"IMI//-'/i(i>n)  , 

2.  H 

and  in  combination  with  (2.S)  this  proves  (2.4). 

S<>  f<ir  wo  have  considered  one  inhomogeneity.  In  the  case  of  A  inhomogeneities  the 
result  follows  from  the  previous  proof  and  a  localization  argument.  □ 

Lemma  2.1  asserts  that  ||A  —  t/,||//i(u)  —  0(e"/2);  from  the  trace  theorem  it  therefore 
follows  that  || U  —  u,  ||//i/J(a$2)  =  0(f"'2).  This  however  is  not  the  best  possible  estimate, 
in  fact  we  have: 

Lk.MMA  2.2.  There  exists  a  constant  C  such  that 

IK  -  "<  ll//1/?(/H2)  <  C'f"{  1  +  ||,2’|l//-i/2(ai2)}  • 

Proof.  Let  U'  and  ir,  be  solutions  to  the  same  minimization  problems  as  L’  and  u , . 
just  with  !,'•  replaced  by  i.\  From  the  Schwarz  inequality  and  Lemma  2.1  we  get  that 

c-LO)  |  I  *,V(ir-.r1)-v(t'-«,),/.l  js  ^,f,l(IK’H//-i/2(,-Jjt)  +  IK'H/z-i/itan))  • 

u 

Since  ir,  £  1’,.  it  follows  from  (2.2)  that 
(  2. 10 )  [  y  V ic,  •  V(  A  —  u, )  <ir  —  0  . 


Integration  by  parts  yields 


(-11) 


I  yVU'  •  V( V  -  U<)  dx  =  J (U  -  u,)0  (Is  . 

u  Oil 


ow 


since  V(',VU')  =  0  in  9  and  7——  =  4'  on  09.  A  combination  of  (2.9),  (2.10)  and  (2.11) 

Di' 


Hives 


|  I (U  -  »<)'.•  dr  |<  Cf"(||v||//-i/2(an)  +  Ill’ll //- ^  ’ 


mi 

and  thus  l>v  taking  t lie  maximum  over  (/'■  subject  to  || C|| //- M'i(mi)  —  1  an<^  /  t’ds  =  0.  we 

on 

H,et  the  assertion  of  tin'  leunua.  Q 

Him  AUK  2.1.  The  estimate  in  Lemma  2.2  follows  directly  from  the  representation 
formula  for  u ,  which  we  shall  derive  in  the  next  section.  Indeed,  it  follows  that  U  -  u(  is 
of  order  f",  near  09,  in  the  +  norm,  0  </?'</?;  however  the  present  proof  is  simpler 
and  we  have  included  it  for  completeness. 

fj3  A  representation  formula.  We  now  proceed  to  find  an  asymptotic  expression  for 
a,.  This  expression  will  involve  tin1  functions  T,  ( j  =  l,...,n)  which  solve  the  exterior 
problems 

AT,  =  0  in  R n\D  . 

T,  =  —ijj  on  OD  , 

T,  is  uniformly  bounded  in  R"\Z?  and 
if  n  >  3  :  T ,({/)  — >  0  as  |  y  |-+  00  ; 

such  <Vj  exist  and  are  unique,  [15].  For  n  >  3  we  shall  also  need  the  function  6,  which 
sat  isfies 


Cbl) 


A 6  —  0  in  R n\D  , 


( -i.2) 


0=  1 


on 


OD  , 


4>(y)  — ►  0  as  |//|  — >  00  , 


The  following  lemma  is  well  known,  and  is  only  stated  here  for  the  convenience  of  the 
leader. 

Ll  M M  A  3.1.  Let  <1»(  // )  he  a  uniformly  hounded  harmonic  function  in  R "\D.  If  n  >  3, 
assume  that  <!•(//)  — >  0  as  |  //  |— >  oc.  Then  there  exists  some  constant  c  such  that 

<1  *(//)  -  r  |  ij  +0(|  //  |1-").  and 
V«K//)  =  V(r  |  ij  I2'")  +  ()( |  y  T")  as  |  y  j->  =c. 

S 


& 


>. 

V. 


w 


3HI 

v; 

1 


1 

1 


r  -.1 

$ 

* 

-V-: 


'■A 

0# 


£ 


Briefly  the  lemma  follows  by  considering  the  Kelvin  transform 


-'(//')  - 


/i  i-  n 


7-777).  .'/  urai 


0 


and  noting  that  the  singularity  of  ‘I>'  at  y'  =  0  is  removable  since' 

1 

,  ,  I  o  (log  - 

W)  - 


O  (log  — )  if  71  =  2 

\y\ 


o  (\y'\2  ")  if  n  >  3. 


From  Lemma  3.1  it  follows  that 


/ 


d$, 

( <I> ,  -  $  j  -rr—  )  ds  — >  0  as  if  -►  oo  , 

Ou  Ou 


1Ij»I=«} 

and  (Jreen’s  formula  therefore'  vields 


(3.3) 


Ou  d$  =  J  ^  ds~ 

01)  OB 


Similarly, 

(3.4) 


/ --I 


0$, 


Ou 


ds,  n  >  3 


OB 


OB 


on 


Notice,  by  the’  maximum  principle',  that  0  <  <^>  <  1  in  Rn\Z?  and 


(3.: 


00 

—  >  0  on  OB 
Ou 


(u  is  the  outwarel  normal  relative  to  B). 

We  ne>w  pre>ve  twe>  crucial  estimates,  needed  for  establishing  the  representation  formula. 
LKMMA  3.2.  For  m lv  v  €  H~1^2(OQ)  there  exists  a  constant  C  such  that 


(3.G) 


Ou, 


Ou 


ds  <  Ccn~l  ,  and 


(3.7) 


Ou, 

'!~0^{X 


-At)  (l*r  =  ((PkY'~)(~k)AVU(:k)  +  o{ e"; 


■'{  /») 


L, 


*51 


/or  1  <  k  <  A,  uniformly  in  {;:*}  and  {pk}-  ^  =  («ij)  is  rhe  STnnm'frir  positive  df'finit 
matrix  with  entries 


"i  j  i  n  i  <s 


/*.&* 

on 


* . 


where  <'  is  the  outward  norma]  relative  to  D  and  6n  =  0  if  n  =  2,  <*>„  =  1  if  n  >  3. 

Proof.  Tin*  fact  that  .4  is  symmetric  follows  from  (3.3).  To  verify  that  .4  is  positiv 
definite  we  compute 


=|fl|l  ft'-  / 

on 

+  6"(]^,hr'U  T, 


with  \  =  and  note  that 


-  f  \  <1*  —  j  |  V\  | 2  dy  >  0  (  n  =  2),  and 

on  R2\n 

*/'&*  +  «/ 

on  an  on 

=  J  I  V\  |2  dy-(  J  V\ Vo  dy)~  (  I 


tin  ough  integration  by  parts  and  use  of  Lemma  3.1. 


Vo  |2  dy)~'  >  0  {n  >  3). 


In  prove  (3.G)  and  (3.7)  consider  first  the  case  in  which  there  is  only  one  mhomogeneitv 
with  :  -  0  and  p  =  1  :  u,\  =  e B .  Let  d  be  a  fixed  positive  number  (chosen  sufficient  1\ 

small).  Without  loss  of  generalitv  we  mav  assume  that 


C  { |.r  |  <  -}  C  { |.r |  <  2d}  C  0  . 


Let  \ ,  be  t  he  s<  >lut ion  to 


V  ■  ( ")  V\ , )  =  0  in  ( |.r |  <  d}\eD  , 

\-  -  1  on  0{(D)  and  \f  =  0  on  ( hr  I  =  d}. 


c(  =  xt,\,n  (n.  constant)  , 


and  introduce  the  functions 


(3.10) 


r.(U  =  -  \,(  r)(r.  -  £•'(())),  r  £  {|.r|  <  </)\ fO 


»: ,e{!l,|  <  i)  \  n. 

f  f 


On  c) D: 


Wt{y)  = 


tr(0)  -  tf(«y) 


=  -  W(0)  •  y  -  |  X)  DiiU(v)yiVi 


(3.11) 


-V(7(0)  •  y  as  e  — ►  0  , 


uniformly  in  C,,  +  ^  ,  for  any  0  <  (3'  <  0.  On  \y\  —  d/e  : 


H-',(y)  = 


u((ey)-U(ey)  ut(x)-U(  x) 


with  .r  =  e /y,  |.r|  =  </.  By  Lemma  2.1 


/  (|V(5^iii[U|^P)*<CV 


(f  <kl<2rf) 


V  ■  (7V— - — )  =  0  ,  in  n\{|x|<d/2), 

e 


we  can  apply  elliptic  estimates  to  deduce  that 


(3.12) 


U({x)  -  U{x)  "=J> 

max - <  Ge  2 

Ul=<t  e 


Therefore 


(3.13) 


Wf(y)  |<  C(  2‘  on  |  y  |=  -  , 


and  at  the  same  time  IT,  solves 


(3.1-1) 


V  ■  ('/(f.v)VU'dy))  —  0  in  {|  y  |<  j}\B. 


K 


The  maximum  principle  shows  that  d\Jdu  has  constant  sign  on  d(eB),  he.,  (3.19)  may 


he  written 


(3.20) 


o  -  L'(0)  I 


ds  — >  0  as  e  — »  0  . 


In  conjunction  with  (3.17)  this  gives 


dll(  a  <?  r 
-a —  ds  <C  , 

av 


which  is  exactly  the  assertion  (3.G). 

In  terms  of  the  solutions  <&j  of  (3.1), 


and  thus 


(3.21) 


Wo  =  £^(0)*,  , 


=  ™dD- 


A  combination  of  (3.17),  (3.20)  and  (3.21)  yields 


(3.22) 


where  A  —  (n,.j  is  the  matrix  with 


r,  d$}  .  , 

<iij=  /  (vjxyi  +  -Q—  yi)  dsy 
an 

=1  B\6„-  j 


and  (3.7),  (3.8)  follows. 

We  now  turn  to  the  case  n  >  3.  The  function 


satisfies 


My)  =  ”J(r)  =  r  2  ,  V  =  |  y 


1  "  f 

V  •  (7(e?y)Vm(?/))  <  7(e?/)(m"  +  - - m')  -  (||V7j|z~)rf  — 

r  r 

d 

<0  if  |  y  <  -  , 

e 


provided  <1  is  sufficiently  small.  Using  (3.11),  (3.12)  and  (3.14),  wc  may  compare  W,  with 
Cm(i/)  for  some  positive  constant  C  (independent  of  e)  and  conclude  that 


(3.23) 


I  ll’.(y)  |<  c  j  y  l'2"  for  y  £  {|  y  |<  -  }\B  . 


d 


Therefore  If,  converges  uniformly  on  compact  sets  towards  H'o ,  the  uniciue  solution  of 


(3.24) 


AH'„  =  0  in  R n\B  , 

Ho  =  -VU(0)  •  y  on  dB  , 
TT'o (y)  ->  0  as  |  y  |->  oo  . 


The  statements  (3.15),  (3.1G)  and  (3.17)  remain  valid  and,  as  before, 


U(  0) 


-d-i 


/ 


7 


dx< 

dv 


ds  — *•  —7(0) 


(3.25) 


d(<B)  dB 

04>j 


J  m0(y)ds 


=  ~7(0) 


/ 

OB 


dv 


ds ; 


only  now  the  right-hand  side  is  in  general  ^  0.  (If  =  Cj  |  y  |2  ,l  +  O  (|  y  |*  n)  as 

|  y  |— >  00. n.  >  3,  then  f  [d't>J)/dv  —  0  if  and  only  if  c;  =  0).  The  function  d\Jdv  has 

i>B 

a  fixed  sign  on  0(tB),  and  by  combining  (3.25)  with  (3.17)  we  deduce  that  the  assertion 
(3.G)  holds. 


From  (3.17)  we  obtain  a  relation  similar  to  (3.22)  with  the  additional  term 

c,-U{  0) 

OUB) 


J  =  lim 

<—o  e ' 


/  l{x)lt  xdSl 


appearing  on  the  right-hand  side.  Using  the  barrier  Cm(y)  we  see  that  the  functions  \e(ey) 
are  all  majorized  by  Cin(y)  on  {|  y  |  <  —  }\B.  Thus  \'((ey)  — >  4>(y)  uniformly  with  its 
first  derivatives  in  {|  y  |<  R}\B1  for  any  R  >  0.  Combining  this  fact  with  (3.25),  we  get 


r  ( 0 )( 


dB  dB  dB 


The  formula  (3.S)  follows  by  use  of  (3.4)  in  the  above  expression  for  .7. 

We  have  thus  completed  the  proof  of  Lemma  3.2  in  c.ase  of  one  inhomogeneitv,  with 
0.  />  1 :  a  simple  rescaling  and  a  translation  gives  the  result  for  the  general  case  of  one 
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inhomogeneity.  If  there  are  more  than  one  inhomogeneity,  then  we  can  apply  the  previous 
argument  to  the  inhomogeneities  one  at  a  time,  provided  we  verify  that 

du( 


0 


( 3.2G )  J  7  ds  = 

Zk  +(pk  13 

for  each  1  <  k  <  K.  To  prove  (3. 20)  we  return  to  the  uf  defined  in  the  introduction;  it 
is  clear  that  (3.20)  holds  for  all  uf,  and  since  ue  =  lim  uf  (in  Hl(Q))  the  same  identity 

G  — *  OO 

holds  for  uf.  [) 

Let  N(-,y)  denote  the  Neumann  function  in  Q  corresponding  to  7,  i.e. ,  the  solution  of 


-V  •  (7VjV)  =  6y  in  n , 

0N  1  00 


dv  I  dQ 


with 


J  Nds  =  0  ; 


an 


i'  denotes  the  outward  normal  to  dQ.  If  y  €  Q\lo(  then 

3it{y)  =  -  J  ««V  •  (7  V  Ar)  dx  =  ~  J 


ds 


;  3.27) 


an 


f  0u(  f  du( 

-  /  N  ds  -j-  /  7^—  N  ds. 


dv 


d*>< 


an 


Since  7( dN /dv )  is  constant  on  dQ  and  f  u(  ds  =  0,  it  follows  that  f  ut~f(dN  /  dv)  ds  =  0. 

an  an 

Combining  this  and  the  identity  ~/(du(/dv)  ~  if  on  dQ  with  (3.27),  we  get 


( 3.2S ) 


«t(l/) 


=  —  J  7  N  ds  +  J  if  Nds  . 


duj( 


dQ 


Note  that  v  denotes  the  outward  normal  relative  to  u>(,  The  first  term  on  the  right-hand 
side  of  (3.28)  may  he  written  as 


K 


(3.29) 


£  / 


*  *a(  +  (pk  13) 


dv 


!< 


“  r  f)  r 

=  -^2N(z,,y)  /  7  ds  —  V  :N(zk,  y)  ■  j 

r  =  i  J  l/  J 

-  / 


dll(  (  \  1 

(x  -  zk)ds j 


+  0(f 


<}{~k  +(pk  13) 
,du 


d(  :k  +cpk  D ) 


n  dx 


ds)  . 


rl(  :k+f  I’k  13) 


15 


\Ye  already  know  that 


■  *(;,-  t'W  ID 

Using  Lemma  3.2  and  (3.28),  (3.20)  we  therefore  get: 


Lkm.MA  3.3.  Their  holds: 


(3.30) 


i< 

v.(y)  =  -fn  T(zk,y).AVU{zk) 

=  i 

+  J  iHx)X(.r,y)  dsx  +  enr/(e,y,  {/>*}. {'*.  })  , 
as  i 


where  the  matrix  A  =  (a,j)  is  given  by  (3.8)  and  //(e,  y,  {/9jt},  {-2*:} )  as  well  as  Vyr/(e,  y,  {^A } ) 
converges  to  zero,  as  (  approaches  zero,  uniformly  with  respect  to  y  £  Cl  ,  {/>*.}  and  {•:/,}. 
provid'd  <!ist  (//,  {~*})  >  d  >  0. 


Note  that. 

J  4'(x)N{x,y)  dsx  =  U(y)  . 

(Hi 

In  general  the  matrix  A  cannot  be  computed  explicitly.  One  exception  is  the  case  when 
B  is  the  unit  ball.  In  that  case  —  ~V)t  I  V  I”  and  (for  n  ^  3)  <t>  =[  y  |2~".  A  simple 
computation  gives  A  =  u where  /  is  the  identity  matrix  and  u>„  —  meas  {|  y  |  =  1). 

§4  Continuous  dependence.  The  Neumann  function,  introduced  earlier  is  symmet¬ 
ric  in  its  arguments  in 

(12  x  £2)  \  diag  (12  x  12).  It  furthermore  has  the  form 


N(z,y)  =  N(ij,z)  =  - 


1 


(1.1) 


2iry(z) 

1 


log  \  y-z\  +B(y,z) 


(n  -  2)«;„7(c) 


y-z\2~n  +R(y,=) 


where  fl{  y,  :)  solves 

-Vy  ■  ( 1  (/y)Vy/?( y,z))  = 


i  (y)-(y-z) 


•1.2) 


-'nt(-)  |  2/  -  -  I" 


,  \  0  t?i  \  1  ,  n(y)  (y-=)-i'v 

l(ll)g —  B(y.z)  -  -T-— —  + - 


On 


|  512  |  w'„*/(;)  I  i/  —  ~  I" 


if  n  =  2 
if  n  >  3  . 

yen, 

y  e  on ; 


II  is  not  symmetric  in  its  arguments.  Since  the  function 


1) 


v t ( // )  ■  (y  -  =) 
\y-  -I" 


1G 


is  in  L'’(ft)  for  any  1  <  p  <  it  follows  from  (4.2)  and  elliptic  estimates  that  /?( - .  ~ )  is 
in  lF2,,’(ft)  for  any  z  G  ft,  1  <  p  <  -j^y.  From  (4.1)  we  thus  get 

N(z,-)  G  VFlp(ft)  for  l<p<-^-,zGft, 

(4.3)  n-  1 

AAA  i  lF2Aft\{A)  for  zeQ. 

It  is  easy  to  verify  that  the  function  R(ij,z)  is  differentiable  with  respect  to  2  G  ft.  In 
fact,  for  any  fixed  vector  a  ±  0,  the  function  RQ(-,z )  =  a  ■  S7:R(-,z)  is  the  solution  to 


(4. 4«) 


(4.41.) 


-  V,(7 (y)VyRa(y,z))  =  -  —  a  •  Vz  (V7(y|'(y  ^ 

wn  V  7(~)  I  y-z  !" 

*.(»-->- -  «-v.r 


7(~)  I  y  -  ~  I" 


y  6  <9ft  , 


with  the  normalization 


/  n„(y,z)  <U,  =  4-  J a  -  Vj  *, 

an  an 

~  ~  (n  —  2)u>„  /a'Vz(^(I) 


) - ifn  =  2' 


if  n  >  3 


The  right  hand  side  of  (4.4a) 

-  £  "  ■  ) - A  "‘AfA1)  ^ 


_J_a  v  •  (y  ~  _  1 

“-’.I  a  A  7(~)  I  V  ~  ~  I"  /  (n-2V„ 

is  in  (ll  1,7(ft))'  ,  n  <  q  ,  since  both 


A^'^'NioAh-'-O  (n- 


log  I  y 


(n  —  2)  and  V;— — - — —  (n  >  3) 

j(z)  \  y-  z\n  2 


are  in  ZTfft),  1  <  p  <  - - -,  and  are  continuous  on  <9ft  (z  G  ft).  The  right  hand  side 

of  (4.41.)  is  continuous  on  dft  (r  G  ft).  By  elliptic  regularity  and  duality  it  now  follows 
that 


Tin'  functions 


Ra(-,z)  G  IF1  ''’(ft)  for  1  <  p  < 

n  -  1 


v.dog  I, ,/-.-!)  =  If-  V'l  ■"  (n=2) 

I  y  -  ~  h 


*V  w.  v.  ,vyvy -yv  WA  »' 


a.VA\y-  =  r")  =  (2-  n)- 


(n  >  3) 


U  ~ 


arc  in  U‘(il)  Imt  not.  in  \V 1  *( ^\{-} )  ■  Consequently,  it  follows  from  (4.1)  and  (4.G)  that 


o 


(1.7) 


AT(  •)  6  L>\il)  for  1  <p< 


n 


n  -  1 

o  •  V;Ar(_~,-)  i  ir1J(Q\{c})  for  :  G  fi 


,  s  6  VL  ,  and 


for  any  vector  o  ^  0. 

Proceeding  once  more  by  the  method  used  to  prove  (4.G)  one  can  show  that 


(4.S) 


D'iR(-.z)  e  Ll’(Q)  for  1  <p<  — 

77  —  1 


Differentiation  of  log  |  (/  —  z  |  give's 


(  D'i  log  |  y-z 


cv  •  /i 

|y--’l2 


2((~  -  y )  ■  a){{z  -  y)  •  ft) 


from  which  we  conclude  that,  in  the  case  u  =  2, 


(D:  log  |  y  -  c  |  a  )  • ;?  is  not  in  Z,1  ( Q\ { r } ) 

for  any  vectors  a  0,  A  0  .  The  same  statement  holds  in  dimension  n  >  3  for  the 
function  |  y  -  z  | 2  —  .  A  combination  of  this  with  (4.1),  (4.8)  yields 

(4.9)  (£>jA’(~,  •)<*)  •  ft  is  not  in  Z1(Q\{r)) 


for  any  vectors  o  0,  A  7^  0  . 
Introduce  the  function  F: 


(4.10)  F({/>k}[\  { ~ it} )(y)  =  k)VzN(zk,y)  ■  AVU{zk)  . 

k=  1 

I  he  Frechet  derivative  of  F  with  respect  to  {/>*}[',  { c is  the  linear  expression 

DF({APk}[\{  Ac,  }{')(£/) 

Is 

-  X  K'(^)(^:A'( rt^MVn-*))  ' 

(111)  *  =  ' 

f  V~,(r,)  •  Ac*)(V.JV(r*,y)  •  A±U(zk)) 
t  ,»J!  ( -  *  )V;  A’(  :k.y)-(  AD]U( zk)Azk ) 

)  :.X(zk.y)  ■  AVU(zk)]  ,  ])  G 


IS 


:*\,*\.,*i 

.v 


Lkmma  4.1.  If 

y  { I  A  i>k.  |  4-  |  A'*.  | )  0  ,  />*■  >  0  and  Vf;(;n)  ^  0  , 

then  the  function  j DF({A/q}{'  ,  {Ac*}{'  )(•)  is  not  identically  zero  in  i  2  \  { c  ^  }  (v  . 

Proof.  If  at.  least  one  Ar*.  is  not  zero,  then  it  follows  from  (4.9),  (4.11)  and  (4.7)  (with 
p  =  1 )  that  DF  is  not  in  Ll (Q\{:k}  {'  )  and  thus  not  identically  zero;  here  we  used  the  fact, 
that  .4  is  positive  definite  and  Vf7(cq.)  ^  0.  If  all  the  A Zk  equal  zero  then  at  least,  one  A pk 
is  not.  equal  to  zero,  and  DF  is  not  in  IF1,1  (fl\{z*}  [' ),  by  (4.7);  thus  again  DF  ^  0.  Q 

For  the  following  lemma  we  need  all  the  assumptions  made  earlier  in  the  introduction, 
in  particular  (1.3), (1.4)  and  (1.15).  Let 

;  put".  Pi  If') 

=I|C(  {«}■{-*))  -  F(f/>i>.(-i))lli»(n 
+  II I;  |F({W). {-•>!)-  F(Pi>.{4>]lk~<n  . 
where  T  is  a  fix(.*d  nonempty,  open  subset  of  Of l. 

Lk.MMA  4.2.  There  exists  a  positive  constant  <5  such  that  if 

then 

(i)  A  =  K'  and ,  after  appropriate  reordering, 

oo  i  i  + 1  />.--/>;  i  <  cH({ptyy  t-nt.  w\\  put  *  ■ 

The  constants  fi  and  C  depend  on  the  same  parameters  as  2>o.f0  and  C  in  Theorem  1.1. 

Proof.  Suppose  the  assertion  K  —  K'  is  not  true.  Then  there  exist  sequences 

_  /,(">  )iA'  „(m>  _ 

—  / i  ’  P  —  \Pk  I i  > 

=  {-T’  'if  ■£,my  = 

with  A  7^  A '  such  that  the  corresponding  H  converges  to  0  as  m  — »  oo  (since  K.K1  < 
M o  =|  f2  |  /  |  {|  ,r  |  <  <7()/2}  |,  we  can  choose  K  and  IF  to  be  independent  of  m).  Passing 
to  a  subsequence  we  get 


.(»») 


.(>»)  ' 


z'  .  p{’"]  -»  P  and  plm)  '  -+  p 
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Shut  //  approaches  zero, 

F(,ki)(u)=  lim  F(p(,"\r( '"’)(//) 

—  tn  —  x 

-  lim  F(/.(",)\r(,">  ')(»/)  =  F(p',c')(i/),  y  €  I\ 

»n  —  x  — 

and  similarlv  for  -7-  F 

-  f#*/ 

O 

t-  F(p,c)(y)  =  —  F{p',:')(y)  ,  y  G  T  . 
cm  -  yi' 

1  lilis  the  function  C(i/)  ■==  F(p,r)(y)  —  F(p',c')(;/)  is  a  solution  of 

v  •  (■>(?/ )VG’)  =  0  in  n\({^}|' U{4}|v") 

which  vanishes  together  with  its  first  derivatives  on  T.  By  uniqueness  of  the  solution  to 
the  Cauchy  problem  ([3],  [7])  it  follows  that  G  is  identically  zero.  i.e.. 

(•1.12)  F(,k  =  )  =  F(p',c')  in  Q  \({:k}['U 

The  components  of  />  and  />'  are  all  positive,  so  F(p,r)  has  a  singularity  at  precisely  each 
of  the  points  c*.  and  F(//.  F)  has  a  singularity  at  precisely  each  of  the  points  :'k  (by  (4.7)). 
Therefore  (1.12)  is  a  contradiction,  and  we  conclude  that  K  =  K' .  The  same  argument 
shows  that,  after  reordering. 

(1.13)  |  —  :[  |  +  |  p*  —  /q.  |— >  0  as  //—>(). 


The  proof  of  (ii)  proceeds  along  the  same  lines.  Assume  that  the  assertion  is  not  true. 
1  hen  there  exist  sequences  r ( C"1*  and  //m)  (each  element  a  A"  — vector)  such 
that 


(  1.14) 


,(»>)  ~("0.  „(">)  '  '  • 


//(p(m),~(m);p( 


0 


Li  (I  ~k  ~~k  I  +  \Pk  -  Pk  I) 

Since  the  denominator  in  (4.14)  is  bounded,  it  follows  that  H  —>  0.  and  thus  by  (4.13) 

i< 

/("')  / 1  ('">  ,(>")'  I  ,  I  (m)  ("i)'  a  n 

<<  ~  l  +  lPk  ~Pk  D-^0' 


after  reordering.  Passing  to  a  subsequence  we  may  assume  that 

lime'”0  =  lim  c.*  ”,)  '  =  c  . 
lim  p*'"’  —  lim p(m!  —  p  , 


lim 


(in)  Ini) 


A:  1  . 


(ml  (in)  ' 

,■  !lk  (Ik  \ 

Inn - - — - -  —  A/i  : 

f/(m)  '—k 
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note  that 


Y  (  I  A'k  I  +  I  A/u-  I  )  =  1  . 

1  '/ 

I 

From  (4.14)  and  tho  definition  of  H  it  follows  that  J 

j 

iim'ua, (Anit  )n 

where  /)/•’  denotes  the  Frechet  derivative  of  F  at  (/>.c).  However.  DF{  { Apr} .  {  Ar;- } )(  y ) 
is  clearly  a  solution  of  the  homogeneous  equation 

Vv  ■  h(y)VJ,(JDF)(i/))  =  0  in 

(see  (4.11)).  and  by  uniqueness  of  the  solution  to  the  Cauchy  problem 

DF({ {Act }(' )(y)  =  0  in  . 


1  his  is  a  contradiction  to  Lemma  4.1.  Q 

/boo/  >>t  Theorem  1.1.  If  ||ti,  —  u\ || l^(\  )  <  6neu  then  the  representation  formula  in 
Lemma  11.1}  gives 


II  ^  t  )  —  •f'l  {/»*,  }  {'  •  {4 }  ['  ) || /.»(!')  <  f  "  ll'b  “  tl'(  ||  Lx  (D  +  ?/(f ) 


\  K'  I  \  K' 


<  <*>0  +  v(f)  . 


i  •  k>"'  r  .  ,  ,  , 

and  since  — —  —  — —  on  1  it  also  follows  that 

Or  Or 


ij;  tali' )-  I;  <  -i(m 


here  i/(0  denotes  the  maximum  of  2(  |  >/(c.;/.  {/<*},  {c*})  I  +  I  Vv»/(f.t/.  {/>*}.  {-*•})  |)  (from 
Lemma  3.3)  over  t/fL.  {/>k}  F  [do-  F(>]  and  {;;■}  Q  Li.  subject  to  (1.3).  In  terms  of  //. 


<  <*> 


if  we  take  <so  =  6/2  and  f  <  e,)  so  that  2i/(f)  <  6/2.  Applying  Lemma  4.2  we  conclude  that 
K  --  K'  and 


=k  -  41+  I  in  -  in  l<  CH  <  Cf  "|l»  -  u',  (D  +  i}(f)  . 


where  i/(c)  —  20/(0  — *  0  as  f  — »  0. 
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R  km  A  It  K  4.1.  If  f*  =2  and  if  U  |(m  has  only  finitely  many  relative  extrema  then 
there  an-  at  most  finitely  many  points  in  ft  where  VU  can  vanish  (cf.  [1]).  If  the  points 
c i , .  .  .  ,  c /x  are  restricted  to  stay  away  from  these  critical  points,  then  the  non-degeneracy 
condition  (1.15)  nuiy  he  dropped. 


R  IM  AUK  4.2.  In  the  special  case  of  a  constant  reference  conductivity  7  and  11  =  2, 
the  results  in  this  section  have  some  similarity  to  the  results  about,  the  location  of  poles 
for  meromorphic  functions,  found  in  [  1 G] .  If  T  is  all  of  <9ft  then  Lemma  4.2  r  >llows  directly 
from  the  representation  formula  in  Section  3  and  the  analysis  in  [  1G] .  On  the  other  hand 
if  I'  is  a  proper  subset  of  <9ft  then  this  analysis  would  give  a  Holder  estimate  and  not  the 
stronger  Lipschitz  estimate  of  Lemma  4.2 


Ri;,M,\HK  4.3.  One  could  raise  the  question  of  “best  possible"  choice  for  T  or  or 
what  are  the  shape's  D  that  are  the  easiest  or  most  difficult  to  find.  This  analysis  would 
be  somewhat  in  the  spirit  of  [10],  although  the  results  would  differ  due  to  the  different 
<  1!  >  )ect  1  v< 's. 


!i 5.  Small  inhomogeneities  with  zero  conductivity.  The  results  of  Sections  1-4 
can  be  extended  to  the  case  where  the  inhomogeneities  have  zero  conductivity,  that  is,  to 
the  case  when'  is  replaced  by  For  n  =  2  this  follows  from  the  fact  that  problems 
( 1.5)  and  ( l.G)  are  very  closely  related  by  duality.  The  flux  a(  =  yVu®  in  Q\u)f  .  a,  —  0 
in  vC,  minimizes  the  functional 

\  1 1-'  k  I2  <!'■ 


subject  to  V  •  rr  =  0  in  ft.rr  ■  v  =  (,>  on  Oil  and  (7  =  0  in  jot.  Assuming  that  ft  is  si 
connected  we  get  that 

,  0  o 

n<  =  v< ■  o —  r<)  ' 

OI 2  U'r\ 


wIhtr  r,  is  the  mmnnizer  o 


\  J  5"1  I  Ve  I2  (/,' 


subject  to  Ve  =  0  in  wC,  and  e  =  'L  on  3ft  ( 'L  is  determined  as  a  clockwise  integral 
of  1,’).  The  Neumann  problem  (l.G)  with  inhomogeneities  of  conductivity  0  is  therefore 
equivalent  to  a  Dirichlet  problem  with  infinitely  conducting  inhomogeneities  and  reference 
<  ondiict  ivity  7  “ 1 .  Since  the  method  of  jjjjl  —  4  applies  equally  well  to  the  case  of  the 
Dirichlet  problem,  we  conclude  that  Theorem  1.1  extends  to  problem  ( 1  .Gh  i.e.,  to  11  (  —  11" . 

Consider  now  the  case  n  >  3.  The  results  of  Section  2  can  easily  be  extended  to 
11,  C1.  To  derive  a  representation  formula  for  a,  =  */*’.  we  start  with 


a,  ill)  = 


OX  , 

~l"i  ~X—  -f 

Uu 


I  v.Xils 


1 


( <-f.  (3.27))  whn-f'  i>  is  the  outward  normal  relative  to  u-\.  Suppose  there  is  only  one 

inlioinoficneity,  with  c  =  0.  />  =  1.  We  assume  that  D  is  star-shaped  with  respect  to  the 
origin  and 

c)y 

(0.2)  vp  ±  0  along  OD  (r  =|  y  |)  . 

a  i' 

Let  ‘L ‘ { // )  and  'I'd//)  denote  the  unique  solutions  to 

W(~,(f//)V'I'X'/))  =  0  in  { |  !/  |<  *}\B. 
d'V 

(3.3)  — — L  =  -vt  on  OD  , 


'I'd//)  =  0  on  {|  y  |=  -} 


A'l',  =  0  in  R”\  D  , 

(3.4)  ±=-viondB, 

Du 

d '.(</)  ->  0  as  |  y  |-»  oo  , 

where  v  is  the  outward  normal  relative  to  B.  Notice  that  the  function  m( rj)  =  rh(r)  — 
r  ?  ,  r  --]  //  |  .  satisfies 

—  <  0  on  OD  (bv  (5.2))  . 

Ou 

and  therefore  C'ni(y)  can  be  used  as  a  barrier  for  the  We  conclude  that  d'J  — *  d' ,  in 
C‘  in  any  set  {[  //  |<  D}\D  as  e  — >  0. 


Consider  the  function 


n*,(!/) 


ut{ey)  -  U{ey) 


)'P'(!/)- 


From  energy  estimates  it  follows,  as  in  section  3,  that 
I  Jl*,(//)  I  <C  I  y 


cr  19 


^  +  ”<b  |  y 


on  |  y 


provided  <  <1H b / C ,  0  <  0  <  1/2.  At  the  same  time 


0  uniformlv  on  OD ,  as  f-»0. 


Then  'fori-.  by  comparison. 


H'.(y)  \<i>\y  +  *  in  {|  y  \  <  ~}\D 


fur  f  sufficient  lv  small.  Since  is  arhifrarv.  we  conclude  that 


ii,  .r  — 


l 'U)  +  f  Ur, ( ())'!' ,( j)  +  n'i(e.x)  ,  xfD((D)  , 


where  >/(f,.r)  — +  0  as  f  — »  0.  uniformly  in  x  G  0(eD)  .  Inserting  this  into  (5.1)  we  get 


«.(//) 


i  ,  ow  0X(x,  y)  , 

/  i(x)l(x) — — -  f/sr 

(■/<) 


(/.sx  +  c"t/(c,j/)  +  /  V*(a-)Ar(.r,  7/)  f/.Sj- 


where  i/(c,;/)  — *  0  uniformly  in  y.  y  G  ft,  |  y  |>  </  >  0  ;  n  is  the  outward  normal  relative  to 
(Ii. 

The  first  integral  on  the  right-hand  side  can  be  written  as 


[  y(fc){/(f-)(A  N)(ez,ry)<h: 

i'll  I 

=  ~  at)(«, y) rfs*  +  0( 


6"  +  I), 


e  U(c)  r(0)  +  X) (0)  +  0(c'2)  and 


—  A'(ec,y)d.s.  =  0  . 
0\'x 


1 1  =  f"  "i'(h)frj,  (O)-Vj-j  (0.  y)  I  ZiVjds;  +  0(cn- 


where  is  the  /  —  //f  component  of  the  outward  normal  n 

Th  e  second  integral  on  the  right-hand  side  of  (5.C)  can  be  written  as 


/-  =  f"  Y,‘> ({))fri-i(())Ab-> ((>•.'/)  /  +  0(c’,+  I)  . 


In  summary 


<■'  t  ( !/ )  =  f"i(0)V,A'(0,y)^Vr(0) 


f  rH-r) A>, 


V )  +  enrj(e,y) 


with  .1  =  (<i,j)  given  by 


a,j  =  J  v^Zj  +  ^j(z))ds:  , 


i >  is  the'  outward  normal  relative  to  D ,  and  ?/(e,  y)  — >  0  as  e  — »  0,  uniformly  with  respect 
//,  |  y  |>  (I  >  0. 

Similarly,  for  the  case  of  A  rescaled  inhomogeneities,  we  get  the  representation  formula 

K 

My)  =  e"X;^7(~it)V,Ar(c,,y)  •  AW(M 

(AS)  *=1 

+  y  tl’(x)N(x,y)dsz  +  e'b/(e,y,  {/»fc}{~*})  ! 

an 

the  only  difference  between  this  and  (3.30)  is  the  sign  in  front  of  the  first  term  and  the 
matrix  (u,y),  which  here  is  given  by  (5.7). 

From  the  definition  of  VI'7,  (5.4),  and  integration  by  parts  (using  Lemma  3.1)  it  follows 

that 

=  |i?im2+  j  |V*|2d~ 

Rn  \/J 

W  ith  T  =  consequently  .4  is  symmetric  positive  definite.  In  general  .4  cannot  be 

computed  explicit  I'  ’  ver  in  case  D  is  the  unit  ball 


*M)  = 


1  —  n  1^1" 


and  we  compute  that 


«,)  =  — — rr  ,  tcn  =  meas  {|  -  |=  1}  . 

[71  1 J 


Our  main  result.  Theorem  1.1.,  now  immediately  carries  over  to  the  case  of  it® . 
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